Level spacing statistics in a randomly-inhomogeneous acoustic waveguide 
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Dynamics of a randomly-perturbed quantum system with 3/2-degrees of freedom is considered. 
We introduce a transfer operator being the quantum analogue of the specific Poincare map. This 
map was proposed in (Makarov, Uleysky, J. Phys. A: Math. Gen., 2006) for exploring domains of 
finite-time stability, which survive under random excitation. Our attention is concentrated on level 
spacing distribution of the transfer operator, averaged over ensemble of realizations. The problem of 
sound propagation in an oceanic waveguide is considered as an example. For the acoustic frequency 
of 200 Hz, level spacing distribution undergoes the crossover from Poissonian to Wigner-like statistics 
with increasing time, as it is consistent with classical predictions via the specific Poincare map. For 
the frequency of 600 Hz, the level spacing statistics becomes non-universal due to large number 
of nearly-degenerate levels whose contribution grows with time. Occurrence of nearly-degenerate 
levels is attributed to bifurcations of classical periodic orbits, caused by fast spatial oscillations of 
^ ■ the random perturbation. 
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Manifestations of classical chaos in quantum motion are one of the most intriguing fields in the quantum theory. It 
is recognized that quantum behavior qualitatively depends on whether the classical counterpart is regular or chaotic. 
Regular classical motion typically corresponds to laminar and coherent quantum evolution. In contrast, chaos in the 
I classical limit gives rise to diffusive wavepacket spreading 0, Q , fast decoherence || , and irregular nodal patterns 
of eigenfunctions In quantum spectra, classical chaos infers repulsion of the nearest energy levels and Wigner- 
like statistics of level spacings, as opposed to Poissonian statistics in the integrable case @. Thus, level spacing 
distribution is a good indicator of chaotic manifestations into quantum motion. Hamiltonian systems usually have 
mixed phase space with the coexistence of regular and chaotic domains. In the semiclassical limit, level sequences 
^| . contributed from regular and chaotic domains are independent, and spacing distribution obeys the Berry-Robnik law 
J> ■ @, deviating from it as quantum corrections grow Q- 

If the quantum system is subjected to a random perturbation, it is reasonable to expect purely chaotic diffusive 
behavior due to the absence of stable domains 0- However, if we deal with finite-time evolution, there can be 
phase space domains where Lyapunov stability and coherent motion still survive. They can be found by computing 
distribution of finite-time Lyapunov exponents in phase space, or finding eigenfunctions of the classical transfer 
operator Q. In this Letter we consider the method based on the construction of the specific Poincare map for a 
typical realization of the random force [H, [l(| ■ We describe briefly the main scenarios of global chaos and examine 
how they are reflected on the quantum level. In particular, we introduce the quantum version of the specific Poincare 
map, and study how its spectral statistics depends on the length of the time interval. It is shown that this dependence 
can be non-trivial when classical chaos is caused by cascade of periodic orbit bifurcations. Sound propagation in an 
oceanic waveguide is considered as an example. 

We begin with the classical level. Consider the Hamiltonian of a unit-mass particle 



H = y +U(q) + eV(q)aq,t), (1) 

where p and q are momentum and position, respectively, U (q) is an unperturbed potential, s <C 1, V{q) is some smooth 
function, and £,{q,t) is a zero-mean random function. Assume that £ is differentiable, stationary, and normalized as 
(£) = 1. If we consider a sing le realization of £, we can regard it as a deterministic function with a broad temporal 
spectrum. According to [jl Il0| . one can seek for those sets in phase space that transform to themselves without mixing 
in the course of evolution from t = to to t — to + r. Such sets correspond to trajectories that stable by Lyapunov 
within the interval [to : to + r] . They can be found using the map 

Pi+i=p(to+r\pi,qi), qi+i=q(t +r\pi,qi), (2) 
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where p(to + r\pi,qi) and q(to + r\pi,qi) are the solutions of the Hamiltonian equations at t = to + r with initial 
conditions p(t ) = p.^ q(t ) = qi. It means that values of p and q on the i-th step of mapping become the initial 
conditions for the i + 1-th step. After the time shift t' — t — to, the specific Poincare map ([2]) becomes equivalent to 
the usual Poincare map corresponding to the Hamiltonian 

B = £ + U(q)+eV(q)i(q,t f ) (3) 

where is periodic function being a train of identical pieces of £(t) 

£(q,t' + nT)=aq,t'), < t! < r, (4) 

with n being an integer. Thus, we replace the original Hamiltonian by a periodically-driven one. This replacement is 
justified by the restriction to examine evolution only within the interval [to : to + t\. The main property of the map 
@: each point of a continuous closed invariant set of the map 0) corresponds to an initial point of a trajectory that 
is stable by Lyapunov within the interval [to : to + t] . The map @ provides a sufficient but not necessary criterion of 
stability. 

There can be two different scenarios of chaos onset for the map © , depending on whether £ varies with q slowly or 
fastly. For slow variations, transition to global chaos can be described by means of overlapping of resonances fill ll^. 
It is convenient to transform variables p — q to the action-angle variables I — $ [l2[ • Resonances are determined by 
the condition 

mr = nT(I = J res ), (5) 

where m and n are integers, and T is a period of unperturbed oscillations. Strength of a resonance can be quantified 
by the modulo of the Fourier-amplitude 
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V£e- i( - mi > +nvt ")dtfdt', (6) 



where v — 2tt/t. It is natural to expect that |-ff m ,n| decays, on average, with m and n, therefore, only some finite set 
of dominant resonances ([5]) should be taken into account. The case of r <C T corresponds to high-frequency driving 
when all resonances are weak and can be neglected. Then the Hamiltonian ([3]) can be reduced to a time-independent 
and integrable one using the averaging method [l3| . In the limit r>Ta distance between neighboring resonances ([5]) 
scales as r -1 , and they eventually overlap completely, i. e., motion becomes fully chaotic according to the Chirikov's 
criterion. However, small islands of stability can persist even on timescales large compared with T [ol. [Tol] . 

Fast oscillations of £ with q (and, thus, with {f) lead to a more complicated scenario. In this case the integration 
over i? in ^ can be carried out using the stationary phase techniques. The respective stationary phase condition can 
be easily derived if we substitute £ = Texp(ia) + c. c, where T is a slowly-varying term, into (JBJ). Then neglecting 
the slow dependencies of V and T on we obtain the desired condition da/dd ± m = 0. If this condition is 
mlfiUcd simultaneously with vanishing d 2 a/d-d 2 (the second derivative of the phase in ([6])), certain resonances ([5]) are 
enhanced singularly, as well as the corresponding derivatives \dH m _ n /dI\ and \d 2 H m , n /dI |. These derivatives quantify 
perturbation-induced variations of the frequency, and their enhancement can cause violation of nondegeneracy, i. e. 
leads to nonmonotonic dependence of frequency of oscillations on the action. This phenomenon is closely linked to 
the occurrence of nonstationary wavelike distortions of the potential, caused by the perturbation imposed. These 
distortions can scatter classical trajectories making them chaotic. As one considers how classical dynamics changes 
with increasing e, the nondegeneracy violation in vicinities of resonances ([SJ corresponds to multiplication of fixed 
points due to the cascade of pitchfork and saddle-center bifurcations. As a result, there arise clusters of splitted orbits. 
If the bifurcations appear, global chaos can emerge even for low r. As r increases, the number of orbits undergoing 
bifurcations grows. A more detailed description of this scenario can 

be found in dlHI. 

One of the most important properties of the specific Poincare map is that area of the overall regular domain weakly 
depends on a realization of £, albeit details of a phase space pattern are typically different. This circumstance validates 
the usage of such a quasideterministic map for studying randomly-driven systems. A more comprehensive analysis 
of the specific Poincare map can be found in [l4T |. In [l6j. the specific Poincare map was used in the presence of 
dissipation. 

Following the analogy with the specific Poincare map, one can introduce a quantum transfer operator F which 
describes evolution of the quantum state between time instants to and to + t under a perturbation involving a 
particular realization of the random function £(t) 

F$(q,t +T)=$(q,t ). (7) 
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It is analogous to the Floquet operator corresponding to a classical Hamiltonian (j3)). One can define eigenstates and 
eigenvalues of the operator F by means of the equation 



F<f> m (q,t 



(8) 



where quantity e m plays the role of quasienergy. In the semiclassical limit K — > 0, quantum behavior described by 
the operator F is consistent with the classical map ([2]). So, one should expect some signatures of the transient stable 
domains in the spectrum of F. We can introduce ensemble-averaged level spacing distribution 
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where s is a difference between two neighboring values of quasienergy, P„(s,t) is a level spacing distribution for F 
calculated with n-th realization of £(q,t). p(s,r) can serve as a representative characteristics describing quantum 
evolution within the temporal interval of a length r. Transformation of p(s, r) with increasing r should be attributed 
to the crossover from an initial laminar phase to intermittency, and, eventually, to ergodic diffusive behavior Q. 

The example we represent in this paper doesn't relate to quantum mechanics directly but is equivalent to quantum- 
mechanical problems in the mathematical sense. That is long-range sound propagation in a randomly-inhomogeneous 
acoustic waveguide in an ocean. The waveguide occurs due to the nonmonotonic depth dependence of the sound speed, 
with a minimum at some depth, and maintains sound waves inside water column with weak attenuation. Wavefield 
is described by the standard parabolic equation 
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(10) 



[Ac(z) + Sc(z, r)] 



Here ko — 2tt f /cq is the wavenumber in the reference medium with c = Co, / is signal frequency. The parabolic 
equation formally coincides with the time-dependent Shrodinger equation for a unit-mass quantum particle. In this 
analogy, range r plays the role of the timelike variable, the inverse wavenumber k^ 1 serves as the Planck constant, 
and U(z,r) is a potential consisted of an unperturbed part Ac(z)/co and a small perturbation 5c(z,r)/cQ caused by 
oceanic internal waves. 

Ray trajectories obey the Hamiltonian equations 



(h = dH_ dp = _dH_ ff = ^. + jjf zr \ 
dr dp ' dr dz 2 
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where p = tan \ is the analog to mechanical momentum, y i s a grazing angle of a sound ray. Perturbation becomes a 
significant factor on long ranges, inducing ray chaos [Til [l7j. We use the following expression for 5c: 



Sc(z, r) 
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Here £ is a small constant taken of 0.0014, B = 1 km is a thermocline depth, and p(r) is a random function being a sum 
of 10000 randomly-phased harmonics with horizontal wavenumbers k r distributed in the range from 27r/100 km -1 to 
2-k/I km -1 , with spectral density decaying as k~ 2 . /i is normalized as (p 2 ) — 1. Despite the model (fT2")l is artificial and 
idealized, it maintains the main properties of realistic internal-wave fields, which are responsible for sound scattering, 
such as the presence of vertical oscillations with a depth-dependent wavenumber. Rate of these oscillations depends 
on the parameter u. We consider two cases: u = 5 (moderate oscillations) and u — 20 (fast oscillations). Unperturbed 
sound-speed profile is taken in the form [l4| 



Ac(z) = C 4- (e- 



< z < h, 



(13) 



where h = 4.0 km is the depth of the ocean bottom, a = 0.5 km -1 , rj = 0.6065, co = 1480 m/s. Ray cycle length (the 
analogue of the period T) in the unperturbed waveguide varies from 37 km to 55 km. Fig. [T] represents examples of 
Poincare plots. In the case of u = 5, the Chirikov's scenario prevails. Area of the regular domain slowly decreases 
with increasing r, and small islands of stability survive for the values of r of hundreds kilometers. When u — 20, the 
scenario of bifurcation-induced chaos dominates, and all islands of stability cease completely already for r = 30 km 
(see Fig. md)). 
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FIG. 1: Examples of phase space portraits constructed via the specific Poincare map. The values of parameters are: (a) u = 5, 
t = fOO km, (b) u = 5, r = 500 km, (c) u = 20, r = fO km, (d) u = 20, r = 30 km. 



In order to examine how peculiarities of classical ray dynamics reveal themselves with finite wavelengths, we solve 
numerically the parabolic equation (|10j) with different realizations of fi, and then construct the respective operators 
F in the basis of eigenfunctions of the unperturbed Hamiltonian H . Reflections off the ocean surface z = are 
fairly modeled using the Dirichlet boundary condition $ = 0. Here we are mainly interested in nearly-horizontal 
propagation, so sound absorption at the bottom can be disregarded, and the bottom is treated as the horizontal 
perfectly-reflecting rigid wall with the Neumann boundary condition d<fr/dz — 0. We calculated p(s,r) for r ranging 
from 5 to 500 km. Number of realizations of £ decreased with increasing r from 1000 to 10 due to the limitations of 
computational facilities. For each realization, the number of levels taken into account is 200 for / = 200 Hz, and 600 
for / = 600 Hz. The dependence p(s) obtained with each r is fitted to the Brody distribution 

P B (s) = 03 + l)A fi sf exp(-ApsP +1 ), (14) 

where Ap — \T((j3 + 2)/((3 + 1))]' 3+1 , and T is the Euler gamma function. The cases (3 = and (3 = 1 correspond to 
Poissonian and Wigner distributions, respectively. Figure [2] shows how the parameter (3 corresponding to the best fit 
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FIG. 2: Brody parameter /3 vs r for / = 200 Hz, u = 5 (empty circles), / = 200 Hz, u = 20 (filled circles), and / = 600 Hz, 
u = 20 (empty squares). 

depends on r. For / = 200 Hz, /3(t) reveals the crossover from Poissonian to Wigner-like statistics. The fluctuations 
of j3 indicate on the fact that, according to the Gutzwiller theory [18(, spectrum of the operator F depends on the 
spectrum of periodic orbits of the map ([2]), which varies with t, and these variations infer some weak nonmonotonicity 
in the behavior of /?. The growth is slower for u — 5 due to the presence of long-living regular domains (see Fig.QJb)). 
For u = 20, (3 increases slower than it is expected by the specific Poincare map. This can be regarded as a manifestation 
of dynamical localization For / = 600 Hz and u — 20 the behavior is different: /3 increases rapidly for low r and 

then, for r > 150 km, decreases. Fig. [3] shows that /3 decreases due to the growth of the number of nearly-degenerate 
levels. This phenomenon can be thought of as the consequence of large fluctuations of spectral density, caused by 
periodic orbit bifurcations This effect is not seen with / = 200 Hz because the majority of bifurcated orbits 

appear poorly isolated and don't give individual contributions into spectrum. For / = 600 Hz and u — 5 the influence 
of bifurcations is much less pronounced, therefore, f3 grows until r ~ 400 km, and only then starts decreasing (not 
shown) . 

In summary, we calculated level spacing distribution for the operator which is the quantum analogue of the specific 
Poincare map and describes evolution of a randomly-driven quantum system between two time instants. If global 
classical chaos emerges solely according to the Chirikov's scenario, one expects the transition from Poissonian to 
Wigner-like statistics with increasing the length of the time interval. As an example, we considered guided sound 
propagation in the ocean, when fast depth oscillations of the sound-speed perturbation give rise to cascade of periodic 
orbit bifurcations. As a result, there occur large non-universal spectral fluctuations, and level spacing distribution 
for high acoustic frequencies doesn't exhibit the transition from Poissonian to Wigner-like statistics, despite of the 
onset of ergodic chaos. However, this transition is recovered for smaller acoustic frequencies due to the wave-based 
suppression of bifurcation influence. We suppose that transformation of p(s) with r can efficiently describe quantum 
decoherence in systems subjected to temporal noise, when Chirikov's scenario should dominate. 

This work was supported by the grants from the President of the Russian Federation (project MK-4324.2009.5), 
RFBR-FEB RAS (project 09-05-98608), and the "Dynasty" foundation. Authors are grateful to A.R. Kolovsky, 
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FIG. 3: Averaged level spacing distributions computed with / = 600 Hz, u — 20, r = 120 km (dotted) and r = 500 km (solid). 
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